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Abstract 

We apply the Lee-Suzuki iteration method to calculate the linked- 
folded diagram series for a new Nijmegen local NN potential. We 
obtain an exact starting-energy-independent effective two-body inter- 
action for a multi-shell, no-core, harmonic-oscillator model space. It 
is found that the resulting effective-interaction matrix elements can be 
well approximated by the Brueckner G-matrix elements evaluated at 
starting energies selected in a simple way. These starting energies are 
closely related to the energies of the initial two-particle states in the 
ladder diagrams. The "exact" and approximate effective interactions 
are used to calculate the energy spectrum of 6 Li in order to test the 
utility of the approximate form. 



1 Introduction 



Conventional shell-model (SM) calculations often assume an inert core with 
a few valence nucleons as active particles. The calculation of the two-body 
effective interaction for the valence nucleons can be conveniently divided 
into three steps: (1) Calculate the Brueckner reaction matrix G |IJ] from 
a realistic NN potential; (2) Calculate the two-body Q-box |2| from the 
G-matrix; and (3) Calculate the folded diagrams || from the Q-box. The 
second step provides major difficulties, because one is unable to evaluate the 
core-polarization diagrams to all orders and there is no sign of convergence 
within the lowest few orders p|. This difficulty might be avoided by the 
use of a no-core model space ||, for which all the nucleons in a nucleus are 
treated as active. In such a model space, because there are no hole lines, all 
the core-polarization diagrams are absent and the two-body Q-box reduces 
to the G-matrix. The folded diagrams can be calculated from the Q-box 
by using iteration methods proposed by Kuo and Krenciglowa || or by Lee 
and Suzuki [Q] . Beyond these issues are the largely unexplored questions on 
the role of effective many-body forces in no-core model spaces. 

In previous works [||, pfl, we calculated the low-lying energy spectra for a 
few light nuclei by employing no-core model spaces. We approximated the 
two-body effective interaction by the G-matrix and neglected the folded dia- 
grams. Because of this approximation, our calculations involved the starting 
energy as a parameter, which, in one study Q, was chosen to fit the nuclear 
binding energy. 

In this work, we will calculate the G-matrix using an improved version 
of the Nijmegen potential (Nijmll) [jl(]] and we follow the Lee and Suzuki 
method j7| to sum the two-body folded diagrams to all orders. The resulting 
starting-energy-independent two-body effective interaction is not Hermitian 
but its non-Hermiticity is found to be extremely small. We obtain the 

(21 

Hermitian effective interaction to be used in SM calculation by taking 
the average of the non-Hermitian effective interaction and its conjugate. It 
is pointed out in Ref.[0] that this is an excellent approximation. 
We will also discuss the choice of the two-nucleon Hamiltonian 
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employed in the G-matrix calculation, which determines the intermediate- 
energy spectrum in the ladder diagrams. Currently, due to uncertainties 
in the optimal choices of one-body potentials and methods for treating the 
spurious center-of-mass motion, there is no generally accepted H^ 2 \ It is 
obvious that the two-body effective interactions v& depend on , so that 
one wishes to employ the which best represents the physics of the two- 
nucleon subsystem in the nuclear medium. We will present a choice which 
is physically motivated yet retains simplicity for calculations. 

We will furthermore introduce an approximation scheme which allows us 
to obtain easily, an effective two-body interaction for a no-core model space 
directly from the starting-energy-dependent G-matrix without evaluating 
the folded diagrams. We will demonstrate that the resulting approximate 
form is an improvement over the procedure we introduced in Ref.|{| and 

(2) 

more closely represents the exact theoretical tr ff . 



2 G-matrix and Two-body Effective Interaction 

Assuming that there are only two-body interactions among nucleons in a 
nucleus, the nuclear Hamiltonian for an ^4-nucleon nucleus can be written 
as: 

a I > '/, - /< \i I +T.Vij, (i) 



(I>- t cm)+£ 

\i=l ) i<j 



1<J 

where t is the one-body kinetic energy, Tqm is the center-of-mass (CM) 
kinetic energy of the nucleus, v is the two-body NN potential. 

The expression for the Brueckner reaction matrix G [jij can be generally 
written in the following form: 

G12M = v 12 + v l2 % 7 T7yrVi2, (2) 

where Q is the Pauli projection operator and u is the starting energy. The 
two-nucleon Hamiltonian represents the dynamics of the two-particle 
subsystem in the nuclear medium generated by the remaining (A-2) nucleons 
(also referred to as spectators). It makes sense to optimize the description 
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of this two-particle subsystem in order to minimize the effects of the many- 
body effective interactions which we plan to neglect. Clearly the role of 
the medium on the single-particle states needs to be included and one is, 
therefore, led to introduce one-body potentials into yielding: 

tf( 2 ) = + t 2 ) + V12 + Vi + V 2 = {hi + h 2 ) + U12, (3) 

where h = t + V is the single-particle (SP) Hamiltonian with t the one-body 
kinetic energy and V the mean field generated by the spectators. 

A physically motivated V(r) could be the Hartree-Fock (HF) self-consistent 
mean field or the phenomenological Woods-Saxon (WS) well with appropri- 
ate depth, width and surface thickness. However, to employ either the HF 
field or the WS well is a very computationally demanding project for a re- 
alistic force v\2- For convenience, we assume that the mean field V can be 
approximated by a shifted harmonic-oscillator (HO) potential (see Fig.l), 
namely, 

V(r) ~ Vf& fted (r) = -V + l -mtfr 2 = -V + « HO (r). (4) 

The use of the HO potential in conjunction of the HO basis not only sim- 
plifies the G-matrix calculation but also facilitates the treatment of the 
spurious center-of-mass motion. 

Note that, as far as the low- lying states are concerned, one need not be 
greatly concerned with the obvious fact that, when r goes to infinity, V(r) 
vanishes while u (r) becomes infinite. Actually, in the low- lying states, 
the nucleons in a nucleus remain primarily within the nuclear radius Ra, so 
the shape of V(r) for large r, say, r > 2Ra, plays a less significant role in 
the bound-state spectrum. 

With the approximation stated in Eq.(^), the two-nucleon Hamiltonian 
(H) becomes 

# (2) * (t 1+ t 2 )+v 12 + (uf -V ) + (u2°-Vo) = (hf°+hf o )+v 12 -2V , (5) 
where /i HO = t + u HO is the pure HO SP Hamiltonian. The corresponding 
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G-matrix (0) becomes 



G 12 (u) K « 12 + »i V-ffl + V + , u ) ff " S (6) 

where to' = w+2Vo- In writing the above equation, we have added a constant 
(2Vo) to both the starting energy ui and the two-nucleon Hamiltonian H^ 2 \ 
Obviously this does not change the result for the G-matrix which depends 
only on the difference between u> and H^ 2 \ 

For no-core model spaces, the starting-energy-independent two-body ef- 

(2) 

fective interaction v^J is the ladder diagram series (G-matrix) plus the folded 
diagrams. If we approximate the G-matrix by G ro (lo') as in Eq.(||), is 
written as 

v is - G RO {lu') + (Folded diagrams). (7) 

The folded diagrams can be evaluated by employing the iteration methods 
proposed in Refs.||, 0]. We will use the "vertex-renormalization" procedure 
of Ref.[^]. The input to this method consists of the G(u') and its derivatives 
with respective to u', whose values are taken at an arbitrary (in principle) 
but fixed starting energy a/. 

We use a HO SP basis with M7=18 MeV and a no-core model space 
containing the first 4 major shells (Os, Op, ls-Od and lp-Of). The Q operator 
in Eq.(|6|) is defined to forbid the scattering of the two particles into an 
intermediate state inside the model space (i.e., Q=0). 

For V12, we adopt a new Nijmegen local NN potential (Nijmll) |10| ], 
which was fitted to the world NN scattering data with a nearly optimal 
X 2 per degree of freedom (1.03 per datum). Other potentials (Reid93 and 
AV18) obtained by fitting the same data have a comparable x 2 an d yield 
similar deuteron and triton properties (see Ref.|]13| for more details). The 
derivatives of the G-matrix are calculated numerically through ninth order 
using 11 sets of G(lo') with u/ ranging from -75 MeV to +75 MeV, in steps 

(2) 

of 15 MeV. In Table 1 under column 'viJ ", we list a few diagonal two- 
body matrix elements (TBMEs) of the resulting effective interaction. For 
comparison and further discussion below, we also present, in Table 1, G- 
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Table 1: Selected diagonal TBMEs for the starting-energy-independent two- 

(2) 

body effective interaction uV. These matrix elements should be compared 
with the G-matrix elements listed in the table for three different starting 

(2) 

energies. Values in italics indicate values in good agreement with . Note 
that 2eg° = 3hQ = 54 MeV, 2eg° = 5hQ = 90 MeV, 2egp = 7ftQ = 126 
MeV. 



State 


(2) 


G(u>' = 30) 


GV = 75) 


G(w' = 110) 


"Approx." 


(0 S ? /2 )W=i 


-8.75 


-8.73 


-9.05 


-9.47 


-8.75 




-11.78 


-11.70 


-14.18 


-20.64 


-11.83 




-3.87 


-3.54 


-3.86 


-4.35 


-3.81 




-2.28 


-1.33 


-2.34 


-4.23 


-2.17 


(0^ /2 ) J = ' T=1 


-1.79 


-1.09 


-1.38 


-1.68 


-1.62 


(0^ /2 ) J=1 ' T =° 


-0.78 


0.76 


0.05 


-0.84 


-0.65 



matrix elements at selected values of uj', and at a state-dependent choice of 
J. 

3 Approximation to 

As shown in Eq. ([/]), for a no-core model space, the (^-independent is the 
sum of the ladder diagrams and the folded diagrams, both of which depend 
separately on the starting energy uj (or to'). When the folded diagrams are 
ignored, as was often done in the past, one approximates v e ^ by G(u> ). 
It is obvious that the contribution of the folded diagrams correlates with 
u. Below, we show that a particular choice of uj minimizes the root-mean- 

(2) 

square (rms) contributions of the folded diagrams to the TBMEs of . We 
further show that a state-dependent choice of ui' yields a remarkably good 
approximation to v^. 

In the full theory of the effective Hamiltonian, one has, in principle, in- 
dependence of the mean field V and of uj. However, in practical calculations, 
to minimize the need to calculate higher-order processes, one wishes to make 
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physically sensible choices for these quantities such as discussed above in the 
case of V. 

We now consider arguments that may be presented to suggest a phe- 
nomenological choice for uj that could also simplify the calculation of the 
effective two-body interaction. It is generally accepted that the starting en- 
ergy uj represents the initial energy E% of the two nucleons in the nuclear 
medium. In G-matrix calculations, the two nucleons are treated as two in- 
teracting particles moving in the mean field V. We can think of the energy 
E2 for a two-particle state, which is predominantly \ab)jx (« and b are the 
HO SP orbitals), as given by: 

E 2 = e a + e b + A, (8) 

where e a and e& are eigenenergies of the one-body Hamiltonian (t + V) . The 
quantity A represents the interaction energy and depends implicitly on the 
two-particle state \ab)j : T- 

When the mean field V is approximated by the shifted HO potential as 
we did in Eq.(||), Eq.(||) becomes: 

E 2 c(e™-V ) + (ef°-V )+A, (9) 

where e„° and ef° are the HO SP energies [ef° = (2n, + Zj + §)M) with 
i = a, b\. 

Since the shifted starting energy u/ used in Eq.(||) for the G-matrix is 
related to the original starting energy to through uj' = lo + 2Vq, we have the 
following equation for uj': 

u'~e% + $° + A. (10) 

Note that although E2 and, thus, u are negative for a two-particle state 
bound in the nucleus, uj' is not necessarily negative, because the quantity 
Vo (approximately representing the depth of the potential well generated by 
the spectators) is always positive. This provides a partial justification to 
the choice of uj' made in Ref.Q, where G-matrices at a positive uj' are found 
to yield the approximately correct nuclear binding energies. 



7 



We, therefore, expect the G-matrix for the starting energy u' given by 
Eq.(|i~0|) to be a reasonable approximation to . This is clearly demon- 
strated in Table 1, where the TBMEs of G{ui r ) for three values of uj' are 

(2) 

listed. Prom Table 1, one can see that the listed TBMEs of v e ^ can be 
approximated by those of Cr HO (c</) at a starting energy u/ given by Eq.(|lO|) 
with A ranging from -25 MeV to -15 MeV (numbers in italic). As we men- 
tioned before, A more or less represents the state-dependent contribution to 
the two-nucleon energy from the NN interaction. 

(2) 

If one picks A for each matrix element such that G{uj ) = VL ff , then 
from Table 1, one can also see that A is generally larger in magnitude for 
the J=l, T=0 channel than for the J=0, T=l channel. This is physically 
sensible since the NN interaction for the former channel is more attractive. 
It is also obvious from the table that A is larger in magnitude for lower- 
lying two-nucleon states (e.g., 0s^ 2 ) than for higher- lying ones (e.g., 0py 2 
or 0d^ 2 ), which accounts for the fact that the NN interaction is stronger for 
the former states. 

The observations above lead us to suggest a simple state-dependent 



choice for uj'. That is, we suggest using Eq.(lO) with A taken as a sin- 
gle constant for all states for simplicity. In this way, we hope that G(uj') 

(2) 

will become a good approximation to v^ s . We see that this is, indeed, 
the case by comparing the "exact" results in the first column of Table 1 
with the "approximate" results [i.e. those obtained with uj' of Eq.(|io|) using 
A = — 21MeV] in the last column of Table 1. 

To further illustrate the difference between a state-independent and our 
state-dependent choice of uj' , we define an rms deviation of the matrix ele- 

(2) 

ments of G from those of the starting-energy-independent as 



S V) = {^E [(ab\v^\cd) J>T - (ab\G(uj')\cd) 



J,T 



1/2 

(11) 



where the summation runs over all the 7V=332 two-body matrix elements 
for the SM space consisting of the first three major shells. Note that our full 
model space contains the first four major shells but we are less concerned 
about the matrix elements involving the highest shell, so we omit them from 
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the definition of the rms deviation. 

For a fixed starting energy u>', 5 (to') is plotted in Fig. 2 as a solid curve. 
From the figure, one sees that the smallest rms deviation of about 0.42 MeV 
is obtained when the starting energy u/=75 MeV. 

In Fig. 2 we show as a straight dashed line, the rms deviation when the 
starting energy u' is chosen according to the prescription in Eq.(^) with 
A=-21 MeV. With this prescription of the starting energy, the G-matrix 

(2) 

approximates v^ s rather well since the rms deviation is only 0.13 MeV. 



4 Applications to 6 Li 

(2) 

We now use the "exact" and "approximate" effective interactions tr ff to 
perform SM calculations for 6 Li. The SM effective Hamiltonian is written 
as 

#sm = ( $> - T CM ) + £ v^(ij). (12) 

\i=l / i<j 

The contribution of the center-of-mass spurious motion is removed by adding 
^{H C m — (with A S> 1) to the above Hamiltonian. This is a feature 

available with the oxbash SM code p4|. 

When one compares Eq.([j]) and the above equation, Eq.(12), one sees 

(2) 

that v^Mij) is in the position of Uj,-. Namely, we are replacing the free NN 

potential by the effective two-body interaction. Here we wish to point out 

that the SP potential (u HO ) was used only to determine the intermediate- 
ly 

energy spectrum in calculating v^Mij) from Vij and the SP wavefunctions 
of the basis space. In principle, one expects some contributions from higher- 
order SP insertions. We have not calculated them in the present investiga- 



tions. In Ref.jlS], it has been shown that higher-order SP insertions have a 
negligible effect in large no-core space SM calculations. 

In Table 2, we show the results for the low-lying energy spectrum of 6 Li. 
The calculations are performed in the same model space for which the G- 
matrices and the effective interaction are calculated. But we only allow up 
to 4ftf2 excitations from the lowest-energy configuration [(0s) 4 (0p) 2 ]. With 
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Table 2: The calculated and experimental low-lying energy spectrum for 
6 Li using the "exact" and "approximate" effective interactions as discussed 
in the text. The results obtained using G (lo') at a constant starting 
energy (a/ = O.OMeV and to' = 75.0MeV) are also listed. For the ground 
state, the absolute energy is given. For the excited states, the excitation 
energies are given. All energies are in units of MeV. Since we have not 
included the Coulomb interaction, the experimental ground-state energy 
shown in the table is Coulomb corrected: —31.996 — -Ecouiomb = —33.996, 
where -Ecouiomb = 2.0 MeV is obtained from a HF calculation with the 
Skyrme 3 interaction. 



Jn{T) 


G(u>' = 0) 


G{w' = 75) 


"Approx." 


"Exact" 


Experiment 




-21.497 


-48.386 


-36.854 


-35.655 


-33.996 


3+(0) 


2.481 


2.200 


1.916 


2.054 


2.186 




2.544 


5.246 


5.168 


4.932 


3.563 


2+(0) 


4.955 


6.638 


6.161 


6.306 


4.31 


2ftl) 


5.660 


8.472 


8.406 


8.125 


5.37 


l£(0) 


7.514 


10.057 


9.438 


9.336 


5.65 




11.295 


15.994 


15.719 


15.372 


(N/A) 
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the HO SP basis that we used (HQ=18 MeV), v cf { overbinds the ground state 
by about 1.66 MeV, as shown in the column labelled "Exact". It should be 
pointed out that this result depends quite sensitively on the HO parameter 
TiQ,. Obviously, this is related to the approximation in Eq.(|j), whose quality 
depends on Ti£l. Anyway, we find that when a HO basis with a smaller 
TiVL is used, the resulting two-body effective interaction tends to overbind 
6 Li by an even larger amount. We notice that this is also a feature of the 
results obtained by Poppelier and Brussaard in Ref. ||l6[| (see Fig. 7 in this 
reference), although in that work, the effective interaction has some residual 
dependence on the starting energy. 

(2) 

The calculated excitation energies shown in Table 2 for ir ff J tend to be 
higher than the experimental results but the experimental level sequence is 



more or less reproduced. In Ref.[16|, the excitation energies are even higher. 

It is not clear to us why the effective interactions obtained through the 
Lee-Suzuki procedure from the G-matrices using an HO SP basis with h£l 
smaller than 18 MeV tend to overbind 6 Li. It is quite likely that for light 
nuclei, the approximation made in Eq.(|j) of replacing the mean field gener- 
ated by the spectators by a shifted HO potential with hO, < 18MeV requires 
significant corrections such as effective three-body forces. Further investi- 
gations on this are necessary. 

Shown in Table 2 under column "Approx." are the results of the SM 

(2) 

calculation using the G-matrix (instead of v^J ) calculated at the starting en- 
ergies given by Eq.([H]) with A = — 21MeV. This G-matrix has been demon- 

(2) 

strated in the previous section to be a good approximation to v e ^ when the 
individual matrix elements are compared (see Table 1 and Fig.l). Appar- 

(2) 

ently it is also a good approximation to v c $ when tested by evaluating the 
energy spectrum. 

In contrast, the G-matrix evaluated at any constant starting energy u/ 

(2) 

is not a very good approximation to v^, s . One can see from Table 2 that 
the calculated ground-state binding energy of 6 Li using G-matrix at to' = 
is about 14 MeV smaller than the "exact" result (-21.497 MeV vs -35.655 
MeV). The calculated energy spectrum is also different. Indeed this G- 
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matrix has an rms deviation of about 1 MeV in its TBMEs from those of 

(2) (2) 

ir ff (see Fig. 2) and should not be expected to approximate well. 

When the starting energy to' is restricted to be a constant, the G- matrix 
at to' = 75 best approximates vzJ when the rms deviation 5 (to') in the 
TBMEs [Eq.(|ll])] is used as the criterion (Fig. 2). The SM results using 
G(to' = 75) are also listed in Table 2. The ground-state energy is clearly too 
low compared to the "exact" result. 

5 Conclusions 

We have succeeded in evaluating a starting-energy-independent effective 

(2) 

two-body interaction tr ff ; for a large no-core model space for the new Ni- 
jmegen potential |]10| . Our main conclusion is that v^g can be well-approximated 
by the G-matrix elements evaluated at a starting energy which depends 
on the energy of the initial two-particle state in the ladder diagrams. We 
have seen from Table 1 and Fig.l that for the effective-interaction TBME 

(2) 

(ab\v eft '\cd)j : x compares favorably with the corresponding TBME of G{u ) 
where to' is given by 

a/ = £° + e?° + A. (13) 

For 6 Li with /iSl=18MeV, we found that A = -21MeV is a good choice. 
Notice that when to' is away from the poles of G(to'), a variation of a few 
MeV in A is not significant as the G-matrix element is a very slowly varying 
function of to' in this case. 

We emphasize that the choice of could be very important, as it 
determines the intermediate-energy spectrum in the two-nucleon multiple 
scattering processes (ladder diagrams). Different will lead to different 

(2) 

v^ s . In this work, we have approximated the mean field generated by the 
spectator particles by a shifted HO potential, which seems to give a reason- 
able description for 6 Li, when the HO characteristic parameter 71^2=18 MeV 
is used. 
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Figure Captions 



Fig.l The phenomenological, the shifted HO and the pure HO potential 
wells. The mean field V(r) generated by the (^4 — 2) spectator nucleons 
may be approximated by the phenomenological well, which is further ap- 
proximated by the shifted HO potential for convenience. In the G-matrix 
calculation, we have added a constant shift 2Vo to both the two-nucleon 
Hamiltonian and the starting energy u, so the pure HO potential is 
used along with a shifted starting energy a/ = u + 2Vq. The same amount of 
shift added to and u clearly does not have any effect on the G-matrix, 
because only the difference u — enters the G-matrix equation (]6|). 

Fig. 2 The rms deviation S(uj') of the two-body matrix elements of G at 
a fixed starting energy a/ from those of the starting-energy independent 
two-body effective interaction , as defined by Eq(|TT|) (solid curve). The 
horizontal dashed line is the rms deviation when the starting energy for the 
G-matrix elements is given by the state-dependent choice in Eq.(^) with 
A = -21MeV. 
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